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Control strategies for dissipative preparation of target quantum states, both pure and mixed, and 
subspaces are obtained by switching between a set of available semigroup generators. We show that 
the class of problems of interest can be recast, from a control theoretic perspective, into a switched- 
stabilization problem for linear dynamics. This is attained by a suitable affine transformation of the 
coherence-vector representation. In particular, we propose and compare stabilizing time-based and 
state-based switching rules for entangled state preparation, showing that the latter not only ensure 
faster convergence with respect to non-switching methods, but can designed so that they retain 
robustness with respect to initialization, as long as the target is a pure state or a subspace. 


I. INTRODUCTION 

State preparation tasks are fundamental building 
blocks in control of experimental physical systems and 
quantum information protocols for quantum technolo¬ 
gies. Quantum error correction implies preparation of 
a “code” of states [ill, US, US, HE; measurement-based 
quantum computing starts from a multi-partite system 
prepared in an entangled state [ 23 ] ; state preparation is 
a benchmark for open quantum system simulations Q; 
even a complete quantum computation algorithm, if im¬ 
plemented via engineered dissipation [4l|], can be seen 
as a state-stabilization problem. If the initial state is 
unknown or the initial state is not unitarily equivalent 
to the target, dissipative dynamics or measurements and 
feedback control must be involved in any effective prepa¬ 
ration protocol. Existing techniques for designing such 
dissipative effects include feedback stabilization I3.l5l.l46j. 
either measurement-based [3, HE HE HO, EH or coherent 
HE H3, S3 , open-loop techniques based on Hamiltonian 
control [311 ] . and environment engineering via both open 
and closed loop techniques [lE [32], HE HE HE; 

In this paper, we develop a new set of techniques, 
based on alternating continuous-time dynamics following 
an open-loop switching law. We shall consider piece-wise 
constant dynamical generators in Gorini-Kossakowskii- 
Sudarshan/Lindblad form. These are a particular case of 
generators of (two-parameter) semigroups of completely- 
positve, trace preserving maps. 

The idea of approximating continuous dynamics via 
fast switching of generators is at the basis of geomet¬ 
ric methods for nonlinear control, and it has been suc¬ 
cessfully employed to explore controllability properties of 
quantum semigroup evolutions in [3, [E ;.E EE] ■ 

However, our perspective is different, as the key prob¬ 
lems we want to address concern finding a switching law 
that asymptotically prepares a target state (or subspace), 
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irrespective of the initial state , as well as studying which 
is the best performing strategy among the effective ones. 
Robustness with respect to the initial state is sought af¬ 
ter in typical experimental and information-processing 
tasks, as those recalled above. 

Switched (and, more generally, hybrid) system stabil¬ 
ity is well studied in classical control theory HE HE] and it 
is in general difficult to assess. In particular, even switch¬ 
ing between asymptotically stable linear evolutions could 
lead to unstable global behavior 20!]. When quantum 
(one- or two-parameter) semigroup dynamics are studied 
as linear systems on operator spaces, two specific features 
become apparent: (i) no unstable behavior can emerge, 
since they are contractions with respect to the trace dis¬ 
tance between operators [IJ ; (ii) in order to preserve the 
trace of the state, they always have a zero eigenvalue, and 
hence admit multiple equilibria on operator space. Sta¬ 
bility and convergence to fixed points in Schrodinger’s 
picture is thus better studied in the so-called coherence- 
vector representation, which leads to an affine form of 
the generators. While the classical results for linear sys¬ 
tems cannot be directly applied, we show that, condition¬ 
ally to the existence of a common fixed point, a common 
change of representation can be constructed so that mul¬ 
tiple affine differential equations can be transformed into 
linear ones. By exploiting this transformation, we are 
then able to derive sufficient conditions for the existence 
of effective, deterministic switching laws that stabilize 
the target. These conditions build on existing control- 
theoretic methods, and rely on the existence of a convex 
combination of the available generators that is stabiliz¬ 
ing. 

We propose two open-loop design approaches. The 
first one is called time-based switching , and cycles 
through a set of generators approximating the stabilizing 
convex combination. The other approach is called state- 
based switching , and requires an estimate of the initial 
state to compute the most convenient generator to im¬ 
plement. This can be done either in real time or off-line 
via simulation, as in model-based feedback design [l3j . 
State-based switching are shown to ensure faster conver¬ 
gence to the target. These methods address stabilization 
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for both pure and mixed states in the same way, and they 
are easily extended to subspace stabilization. 

In principle, since the second group of design methods 
requires the initial state to decide which generator to em¬ 
ploy, it is in general not robust with respect to initializa¬ 
tion. Nonetheless, we show that robustness with respect 
to initialization is retained if the target is a pure state or 
subspace, and the initial estimate is full rank. While in 
this case we have no guarantees on a faster convergence, 
numerical simulations comparing the performance of the 
methods indicate that the advantage persists. 

The paper is structured as follows: We introduce and 
formalize our problem in Section El after reviewing ba¬ 
sic definitions and properties of quantum Markov semi¬ 
groups. Section IIIII illustrates our method to jointly 
transform the affine dynamics into linear ones, so that 
classical control methods of time-based switching can be 
applied. The state-based methods are described in Sec¬ 
tion IIVI and relevant applications to entanglement gen¬ 
eration in multi-qubit systems are presented in Section 
m A discussion of the results and potential developments 
concludes our paper. 


II. SWITCHING QUANTUM MARKOV 
DYNAMICS 

A. Preliminaries: Quantum dynamical semigroups 
and their stability properties 

We here recall the notation and some basic facts re¬ 
garding the dynamics of open quantum systems inter¬ 
acting with a Markovian, yet possibly time-dependent, 
environment. Let Us and %£ denote the Hilbert spaces 
associated to the system of interest S and the environ¬ 
ment £, with dim^s) = N < oo. The joint Hilbert 
space Hss is formally described by the tensor product 
Hs ® T-Ls- Let IB {'H) denote the set of linear opera¬ 
tors on a Hilbert space H. A state for a quantum sys¬ 
tem is associated to a density operator, that is, a trace 
one, positive-semidefinite operator p G ®(H) = {p G 
( B('H)\p = pi > 0, trace(p) = 1}. We will say that the 
support of p is Hp if range(p) = TL p . The evolution of 
a density operator p on kiss is unitary, as the relative 
joint system plus environment is closed. Let us assume 
as usual that p(fo) = Ps(io) ® Pf(to), then at each time 
t > to the reduced dynamics on the finite-dimensional 
system of interest (obtained by averaging over the en¬ 
vironment degrees of freedom by means of partial trace 
23), is given by a Completely-Positive, Trace-Preserving 
(CPTP) linear map £t,t 0 - 

Such kind of evolution is, in general, non-Markovian 
and cannot thus be described by simple dynamical equa¬ 
tions. Nevertheless, in various cases of interest (see e.g. 
[1,1, 0), Markovian models can be reasonably used for 
describing a quantum system. A dynamical system for 
which £t 2 ,t 0 = £t 2 M °£ti,to, f° r an y intermediate time 
t\ , with both the factor maps being CPTP, is said to 


be a two-parameter semigroup [H, |48| . It can be proved 
Guam that the generator of the dynamics J0t(p) can 
be put in symmetrized (Lindblad) form: 

j t p{t)=Ct{p) ( 1 ) 

= -i[H,p(t)\ + J2 L k(t)p(t)Ll(t) - ^{L k (tfL k (t),p(t)}, 

k 

where L k (t) G IB('H) are called noise operators and the 
(effective) Hamiltonian H is Hermitian. In general at 
most n 2 — 1 operators are needed for such identifying a 
generator. If the maps £t.t 0 actually depend only on the 
difference t — to, the generator 0 is time invariant, and 
the one-parameter semigroup is simply called a Quantum 
Dynamical Semigroup (QDS) [lj. 

Let us now recall some basic notions related to invari¬ 
ance and stability of QDSs. A set of density operators S 
is invariant if Vp(fo) G S, 

p(t) = £ t ,t o (p(t 0 )) £ S, Vt>f 0 , 

or simply £ tjto (<S) C 5, Vf > to- Recall that the trace 
distance (or quantum total-variation distance) between 
two density operators p and r is: 

D{p,<j) = ^tr[|p-r|], (2) 

where the operator |A| = VaTA denotes the positive 
square root of A^A. By defining the distance of a state 
from a set as 

D(p, S) := inf D(p,r), 

rGo 

we say that a set S is (simply) stable if it is invariant 
and Ve > 0 there exists <5 > 0 such that D(p(to),S) < 6 
implies Vt > <o : 

D (£t,t 0 (p)iS) < £• 

A set S is Globally Asymptotically Stable (GAS) if it is 
simply stable and Vp(fo) we have 

lim D(£ tto (p), S) = 0. 

£—>• + 00 

Notice that having a GAS state is equivalent to say 
that the QDS is mixing in the mathematical physics or 
Markov chains theory. We prefer the dynamical sys¬ 
tem/control theory jargon since it allows to indicate 
asymptotic properties of both single states and sets. 

It is worth recalling two fundamental results regarding 
stability of QDSs: if a QDS admits a unique fixed point, 
then that state is GAS. In the mathematical physics lit¬ 
erature, this is typically rephrased as “primitivity implies 
mixing” [47j] . Furthermore, it is well-known [2^, [2(| that 
any quantum channel is a contraction in trace norm, that 
is 

D (£t,to(.T),£ t ,t 0 (p)) < D(r,p), Vt>t 0 . (3) 

This entails that any invariant set is automatically simply 
stable. Hence, our focus will be on asymptotic stability. 
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B. Switching dynamics and stabilization task 

In order to implement dissipative control, we assume 
that we are able to “manipulate” or engineer the en¬ 
vironment of the system of interest in order to enact 
one among a set of dynamical generators, and switch 
between these. More precisely, we consider evolutions 
driven by equations of the type m where the operators 
are piece-wise constant, i.e. they are allowed to change 
only in a countable and unbounded sequence of switch¬ 
ing times , say ti,t 2 ,t 3 • ■ ■■ Call C(tj) the (constant) QDS 
generator of the form CD active on the switching inter¬ 
val [tj,tj+ 1 ]. Consider a generic time t and let t k be the 
largest switching time such that t > t k - Since for each 
time-homogeneous generator the evolution map can be 
computed through the exponential of the generator, the 
global dynamics can then be written as: 

£t,i 0 = £t,tk ° £tk,tk-i ° ■ ■ ■ ° £ti,to i (4) 

_ e £(U)(i— tk) e C(tk-i)(tk—tk-i) g£(*o)(ii-io) 

We are now ready to state our state stabilization prob¬ 
lem: 

Problem 1 Given a target set S C and a finite 

set of Lindblad dynamics {£j}fLi such that e^S C S 
for all j and all t > 0, we want to find a piecewise- 
constant switching law j(t), 

j : [0, +oo) —> {1,2,..., M}, 

that admits a countable set of switching times 0 ,ti,f 2 , 
i.e. discontinuity points, so that S is made GAS by se¬ 
lecting the corresponding C(t k ) = £j(t k ) on [tk,tk+ 1 )- 

That is, for any initial state r, 

lim D(£ t t = 0, 

t—>oo 

with St t 0 = gCk-^tk—tk-i.)) * * t g^j(t 0 )(^ i— *°) e 

Particular cases of interest are associated to S being 
a single state p , or S being the set of states that have 
support only on a given subspace H 1 C Hs- Of course, 
the main problem will be determining when such a prob¬ 
lem admits a viable solution, depending on the available 
generators. 

Consider the case in which the target consists in a sin¬ 
gle state p. A trivial case is when there is a generator 
for which the target p is the unique invariant state. By 
the results we recalled before, the target is made GAS by 
simply selecting that generator for the whole evolution. 
An analogous argument holds for subspaces [33]. The 
interesting and more challenging case is when all genera¬ 
tors have the target p as a fixed point, but each generator 
also has other fixed points. Clearly, if there were more 
than one common fixed point to all generators, stabiliza¬ 
tion could not be achieved. A sufficient condition for 
asymptotic stabilizability is provided in the next Section 
(Assumption [[]). 


III. JOINT LINEAR REPRESENTATION AND 
TIME-BASED SWITCHING 

A. Coherence-vector representation 

In order to better adapt classical switching systems 
techniques to the quantum framework, it is useful to 
employ their vectorized form. We can univoquely asso¬ 
ciate a TV 2 -dimensional vector v p to a TV x TV-dimensional 
density operator p by choosing a basis of TV x TV ma¬ 
trices. We choose the so-called coherence vector repre¬ 
sentation: we consider the opportunely scaled identity 
matrix F 0 = jjIn and the set of TV 2 — 1 extended Gell- 
Mann matrices {Fj} as an orthogonal basis of the Hilbert 
space of Hermitian matrices flj] . The components of v p 
can be computed as the (Hilbert-Schmidt) inner product 
v p,j = trace(pFj). Thus, the first component of v p is the 
same for each N-dimensional density operator p as the 
trace is constant: 

A) = ^Tr(p/jv) = 1. (5) 

The dynamics are associated to the vectorized generator, 
that is a TV 2 x TV 2 real matrix whose elements are 

{L} jk =TA{C{F j )Fl\ j, fc = 0,1,.. .TV 2 — 1. (6) 

This is sometimes called the “super-operator” form of 
the generator. Considering that the trace of a density 
operator is a constant of motion, the linearized form of 
(P) reads: 



r o 

0 ... 01 


r 1 1 


Vp = Lv p = 

b 

A 


v r 

(7) 


If we consider only the evolution of the coherence vector 
v r , we obtain an affine representation of the dynamics: 

v r = Av r + b. (8) 

B. The reference case: Unital generators 

Let us suppose, for now, that all the generators Cj are 
unital, i.e. they all share the p = -AI as a fixed state. 
In this case, all the affine components bj of ([5D are null 
and the coherence vector dynamics are therefore linear. 
They can be represented as: 

i) r = AjV r . (9) 

Switching linear systems in the classical case have been 
widely studied and a wide spectrum of results are avail¬ 
able. 

Recall that a square matrix A is said to be Hurwitz (or 
stable, or asymptotically stable) if its eigenvalues have all 
strictly negative real part Hzl and marginally stable if 
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they have non-positive real parts, and some real part 
is zero with corresponding multiplicity one. It is then 
easy to show that an autonomous, time invariant linear 
system x = Ax admits x = 0 as a GAS state if and only if 
A is Hurwitz. We introduce next a typical assumption on 
the switching linear systems associated to matrices {A,-} 
that ensures stabilizability, and that will also represent 
the cornerstone of most of our results. 

Assumption 1 There exists a convex combination A c of 
matrices Aj which is Hurwitz, that is 


C. Joint linear representation for fixed-point 
stabilization 

Given a set of marginally stable Lindblad generators 
that all share the same fixed state p, we would like to find 
a switching law that makes such state asymptotically sta¬ 
ble. As the target state can be different from the purely 
mixed one, we must suppose to deal with non-unital gen¬ 
erators Cj, whose dynamics is univocally associated to 
super-operators Lj and thus to an affine reduced repre¬ 
sentation 


3 S.t CZjAj = A c , = 1. (10) 

3 =1 3=1 


and A c is Hurwitz. 

It is proved [H, E3, EH that if the A'-s are the genera¬ 
tors of linear dynamics, this assumption gives a sufficient 
condition for the existence of stabilizing switching rules. 
Let mod (t, e) be the remainder of the division of t by 
e, namely we have t = q x e + mod (f, e) where q is the 
integer part of t/e. 


Proposition 1 Suppose that Assumption [7] holds for 
some {aj}. Then v = 0 can be made GAS for the 
corresponding switching system © for sufficiently short 
switching times. In particular, there exists e > 0 such 
that the switching law 


a(t) = 


'1 

2 


if mod (t, e) £ [0, one) 
if mod (t,e) £ [a\£, a\£ + a^e) 

( 11 ) 


m if mod (t,£) £[(Y}™ = i aj)£,£), 


makes p GAS. 


i v = AjV r + bj, 


where bj yf 0. In this case we cannot directly use the 
known results for linear switching systems. Neverthe¬ 
less, the following proposition allows us to overcome this 
problem. 

Proposition 2 If all the generators Lj share the same 
steady state v = [1 v r \ T , then there exists an invertible 
matrix 



r l 

0 ... 01 

T = 

Tq 

Tr 



-1 


( 12 ) 


that removes the affine component bj in ©, that is, that 
makes all the generator matrices L[f block diagonal: 


L d 3 =T 


0 


0 ... 0 


Aj 


T = 


0 


0 ... 0 


An 


. (13) 


Furthermore, there exists such T that satisfies: 

1 

T ~ V = 0 


(14) 


The results holds for any controlled switching linear 
system that satisfies Assumption 1. The proof is included 
for completeness in Appendix lAl 

If Assumption [I] holds, the stabilization of the com¬ 
pletely mixed state can be easily obtained by properly 
switching between unital generators. However, the order 
of application is irrelevant, as in the first order approx¬ 
imation the effective generator is the weighted sum of 
single ones. The key point is that, on a cycle long e, each 
Aj must act for a time interval A tj = ctj£, where aj is 
the coefficient of Aj in the Hurwitz convex combination. 

The general problem of stabilizing other states, dif¬ 
ferent from the purely mixed one, requires the use of 
non-unital generators. Nevertheless, we will prove that 
this problem can be solved by formally recasting it into 
the special case of unital generators. In particular we 
will need to find a change of basis such to make all the 
generators of the reduced dynamics © linear. 


Proof. Let us consider the corresponding block struc¬ 
ture of an invertible matrix: 



r Ts 

T P 

T = 

Tq 

Tr 





where T s £ R, T P £ [R 1 ^ 2 " 1 ), Tq £ [R^ 2 - 1 )* 1 , T R £ 
R (JV 2 -1)x(JV 2 -1)^ 

From (fl3l) we have 

TLj = L d T. 

Developing the left-hand side of last equation we get 


ITs 

Tp 1 


r o 

0 ... 01 


■ Tpbj 

TpAj 

Tq 

Tr 


b o 

Aj 

= 

Tpbj 

TrA 3 
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while for the right-hand side we obtain: 


T S 

Tp 1 


0 

0 ... 0 

Tq 

Tp 

= 

AjTq 

AjTp 


(17) 

In order for (flbl) and (El to be equal irrespective of 
bj,Aj, it must be: 

T P = 0. (18) 

Choosing for simplicity Tg = 1 and using 1151) . we get: 



r 1 

0 ... 01 

T = 

Tq 

Tp 



J 


(19) 


which is a change of basis matrix if Tp is invertible. By 
exploiting the lower-block-triangular form of T^ 1 , and 
thus of its inverse, it is straightforward to show that 



0 

0 

Lj = TLjT- 1 = 

T R (bj - AjT-'TQ) 

TrA 3 T~ 1 


We need to show that a Tq and an invertible Tp exist, 
such that Tp(bj — AjT^Tq) = 0 independently of j . In 
order to do this, recall that v r is a common invariant 
state, so that AjV r + bj = 0 for all j. If we choose an 
invertible Tp and define Tq = —Tpv r , we have 

bj - AjT^Tq = bj - AjTf^t-TpVr) = bj + AjV r = 0 . 

□ 

This change of basis allows us to transform all the affine 
dynamics in linear ones. It is worth remarking that, while 
T acts linearly on the vectors v p , it acts as an affine 
transformation on their coherence vectors v r . In fact, 
defining v r as the transformed “coherence” part of v, we 
get v r = -^=Tq + Tpv r , and hence: 

*" 1 T^T Q = T^v r + ^=v r . ( 20 ) 

All the transformed generators Aj behave exactly as 
unital generators. It is therefore possible to apply the 
results we derived in the last section and the existence 
of a Hurwitz convex combination of all the {Aj} is again 
a sufficient condition for the existence of a stabilizing 
switching law. By applying the time-based one described 
in dill) of Proposition [lj we can make the shared, general 
fixed state p asymptotically stable. The final results is 
summarized in the following: 


u r — x R Ur 


Vn 


Theorem 1 Suppose to dispose of a set of generators Cj 
that share a common steady state p. Then there always 
exists an affine transfoi'mation that makes all their coher¬ 
ence vector representations simultaneously linear. That 
is, by (l20l) . we get: 

v r = AjV r . 

If Assumption^ holds for the transformed generators Aj, 
then a switching law as in El makes v asymptotically 
stable, for sufficiently short switching intervals. 


D. Joint linear representation for subspace 
stabilization 


In the previous section we found a class of transfor¬ 
mations that allow to recast the switching stabilization 
problem for general dynamical generators in the easier 
unital case, where we can use the linear system approach. 
Another case of interest in which it is possible to recast 
the problem as a reduced linear system stabilization is 
the asymptotic preparation of subspaces. In the nota¬ 
tion we used in Problem[Q we here consider the target set 
S C 0(H) to be a set of density operators that have sup¬ 
port on a certain subspace Us of H. This class of prob¬ 
lems is of key interest for initialization of error correcting 
and noiseless quantum codes [l^ 
to degenerate ground subspaces 


19 

36 


21 , 

37] 


l49j |. cooling 
and multipar¬ 


tite state preparation using quasi-local resources [36|, |36 
Choosing an opportune basis that respects an orthogonal 
partition H = Hs®Hr, any operator can be represented 
as block matrix: 


X = 


X s 

Xp] 

•o 

Xp\ 


€ V(H). 


( 21 ) 


In order to preserve invariance of S, as assumed in our 
Problem, any available generator Cj must satisfy [33]: 


C; 


PS 

O' 


£s,j(ps) 

O' 

0 

0 

) 

0 

0 


( 22 ) 


for some reduced generator Csj and any density 
(sub)matrix p$ on Hs- 

In order to use linear switching techniques for mak¬ 
ing Hs attractive, it is convenient to employ a vec¬ 
torized form that explicitely separates the components 
related to the elements of an operator basis of 93(%) 
that have support on Hs, and the rest. By writing 
® (H) = ®( Hs ) © ( B(Hs) ± and considering a basis that 
respects such orthogonal subdivision, we obtain a block- 
vector representation for p\ 


1 


VS 


(23) 


where the upper two blocks are associated to the invari¬ 
ant states. The super-operator form of the generators 
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must respect the invariance property (1221) . namely the 
Lj , must have the following block-form: 



■ 0 

0 

0 

4 = 

bj 

4 

L \ 


0 

0 

L i\ 


This simultaneous upper-triangular block representation 
of the generators will play the same role of the trans¬ 
formation T in Proposition [2J allowing us to use the 
same sufficient condition that we used in the unital state- 
preparation case. 

Our aim is to find a switching law for which Vg tends 
to zero for t —> +oo, for any initial state. While the dy¬ 
namics of Vg influences that of vs, the converse is not 
true: thanks to the triangular structure of (Ell) , the evo¬ 
lution of Vg is determined only by the Lj- of the currently 
selected evolution. Thus, if we assume that there exists 
a convex combination of { Lj } that satisfies Assumption 
1, namely 


3ai...a m s.t ^^ctjLj- = Ljr, ^jj ctj = 1, (25) 

j'=i 3 =i 


and Lj- is Hurwitz, then a time-based switching law can 
be derived in the same way as for single states stabiliza¬ 
tion. 


Proposition 3 Let’s suppose to dispose of a set of m 
generators with a common invariant subspace and that 
equation © holds. Then there exists a small enough e > 
0 such that the switching law ED renders the common 
invariant set S = {p € 'D('H)\supp(p) C Tig} GAS. 


E. Switching stabilization with Hermitian 
generators 


In all the previous cases we need to find a Hurwitz 
convex combination of given matrices in order to derive 
a stabilizing switching law. Unfortunately, this task has 
been proved to be a NP-hard problem @ in general, for 
which practical, scalable algorithms are unlikely to be 
found. However, the subclass of dynamics of interest 
here is easier to be studied, since they are all stable. 

There is, however, a class of problems for which it is not 
a necessary to find a Hurwitz convex combination of the 
generators. If all the generators Lj are unital, have her¬ 
mitian symmetry and share only one steady state, then 
it is sufficient to employ each generator infinite times, 
for example by periodically repeating a sequence that in¬ 
volves all the generators. If such matrices are marginally 
stable, as we suppose, there exists a change of basis that 
shows that the euclidean norm of the coherence vector 
can never increase, and switching between them leads to 
the shared steady state. This fact is proven in the fol¬ 
lowing proposition. 


Proposition 4 Consider a set of unital symmetric gen¬ 
erators with linear representation Lj, which share only 
one steady state v. Then there exists a switching law 
that makes v GAS. 


Proof. According to the hypotheses, we are dealing with 
generators which have all the following form: 


0 

0 ... o' 

1° 


0 

Aj 


(A?) 1 ', j = !,••• ,m. 


As the shared steady state is unique, it must be the com¬ 
pletely mixed state, whose coherent part is 

v r = [0 ... 0] T . 

Being symmetric, each matrix Aj is orthogonally diago- 
nalizable with a unitary matrix Uj\ 

A d j = UjAjUj, 

where Uj is an orthogonal matrix. Moreover, being 
marginally stable, Aj has real negative or null eigenvalues 
on its diagonal. Consider the positive-definite Lyapunov 
function 

V(v r ) = Ijvjvr. 

The latter is invariant for unitary change of basis, so it 
can also be expressed as 

V(v r ) = Jj)(UjV r ) = ^v r = V(y r ), 

where we now denote by v r (t) = UjV r {t), with j iden¬ 
tifying the active generator driving the evolution. The 
derivative of the Lyapunov function along the systems’s 
trajectory is 

V(v r ) = Vr~^j- = vjAjVr < 0 , 

as Aj is negative semidefinite. If v r belongs to the 
eigenspace relative to an eigenvalue equal to zero, then it 
is left unchanged. However, if e.g. all the generators are 
periodically employed, as in m with any convex com¬ 
bination with non-zero coefficients, for any v r ^ 0 there 
exists at least one generator A~- for which 

V{v r ) = Vr^- = {v r ) T Aji) r < 0, 

for each v r . Indeed, if this was not the case, there would 
be a nonzero state for which 

d'V'f 

—L = AjV r = 0, j = 1,2,... ,m. 

That is impossible as the origin is supposed to be the only 
shared steady state. Thus the norm of the reduced vector 
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v r keeps decreasing and the density operator p = jjIn is 
asymptotically stable for the switching dynamics made 
by the Lindblad generators Cj that correspond to Lj. □ 

It is worth noting that, as for Proposition [I] the result 
holds for general switching of symmetric linear systems. 
It is the symmetry of the generators that ensures real 
spectrum and thus monotone convergence to the origin. 
The assumption of symmetric generators entails however 
“more” robustness than in the previous case, where As¬ 
sumption 1 was invoked: in fact, there is no need to 
consider short switching time intervals. This also implies 
that, for example, one could randomize the activation of 
the generators, instead of cycling periodically. As long 
as all generators appear for a total unbounded subset of 
times, convergence is guaranteed. 

Of course, the same result applies when all the gener¬ 
ators Aj are symmetric after having been transformed in 
a jointly-linear form. In this case any switching law that 
employs all the generators enough times leads to the sta¬ 
bilization of the only shared steady state, which can, in 
general, be different from the purely mixed one. Hence, 
we can state this corollary of the previous proposition. 

Corollary 1 Consider a set of generators with linear 
representation Lj, which share only one steady state v. 
If there exists a jointly-linear, coherence-vector represen¬ 
tation in which Aj = A t for all j, then there exists a 
switching law that makes v GAS. 


IV. STATE-BASED SWITCHING 
A. Preliminaries 

So far we have illustrated switching techniques for sta¬ 
bilizing quantum states that do not require any knowl¬ 
edge of the state of system, either the initial or the cur¬ 
rent ones. Hence, they are intrinsically robust with re¬ 
spect to initialization. Nevertheless, such switching laws 
do not allow for a straightforward optimization of the 
convergence rate to the target. This is quite clear if we 
consider any convex combination of two matrices describ¬ 
ing linear systems: 

A c = ol\A\ + 02^2, ax + a 2 = 1. 

If the switched system is stabilizable, there exists, in gen¬ 
eral, a range of values for a\ (and then for 0 * 2 ) that makes 
A c Hurwitz. One could in principle try to compute the 
spectral abscissa, i.e. the real part of the “slowest” eigen¬ 
value and use that as a worst-case estimate of the asymp¬ 
totic speed of convergence. However, the dependence of 
the spectrum on the on coefficients is in general hard 
to establish, and it is becoming increasingly clear that 
the spectral properties of evolution, also in the quantum 
case, may not be the only interesting performance index 

Hii 


We hereby present a way to circumvent these difficul¬ 
ties that let us design a switching strategy with steepest 
descent, or at least a guaranteed minimal decrease of a 
quadratic function (for quantum systems, a variation of 
the Hilbert-Schmidt distance). In order to do that, we 
need to have an estimation of the initial state of the sys¬ 
tem. We can then exploit this additional information 
in order to find control laws that optimize convergence 
rates. 

The method follows closely the ideas presented in [29| 
for general linear systems, that we revisit in the following. 
As we did before, we introduce the strategy for a generic 
set of switching and marginally stable linear dynamics 
x = AjX. Again, we need to suppose that Assumption [T| 
holds, namely that there exists a Hurwtiz convex com¬ 
bination A c = y . cujAj of the linear generators of the 
vectorized dynamics. 

First, we associate to the asymptotically stable system 
generated by A c a positive-definite Lyapunov function, 

V{x)=x T Px , (26) 

where P is the only positive solution of the Lyapunov 
equation, A (T P + PA C = — I. According to standard lin¬ 
ear Lyapunov theory (see e.g. [l7j|L this ensures that V 
decreases along the trajectories generated by x = A c x 
with rate of — ||:r|| 2 . The function itself can be used as 
a pseudo-distance with respect to the target, and its de¬ 
crease is a good estimate for the convergence rate. 

Next, take this P and, for each marginally stable Aj, 
compute the symmetric, but not strictly negative, matri¬ 
ces Qj = Aj P + PAj. Notice that this allows to easily 
determine the derivative of V along the trajectory driven 
by Aj, since: 

-y- V(x ) = x T (A^ P + PAAx = x r QjX. 

dt J 

Then the following Lemma holds: 

Lemma 1 For each state x 7 ^ 0 there exists a subsystem 
j for which 

x{t) T Qjx{t) < —x(t) T x(t). 

Proof. As A c is a Hurwitz matrix, according to the 
definition (l26l) of P, 

x(t) T (A^P + PA c )x{t) 

= x{t) T ((^ c*jAj)P + P(J2 a jAj))x(t) 
j j 

= axx(t.) T Q\x{t) + ... + a m x(t) T Q m x(t) 

= —x(t) T x{t) < 0 . 

As the sum is negative, then there must exist at least 
one addend which is strictly negative too. Consider the 
j corresponding to the smallest one. Since a convex com¬ 
bination of the rates x{t) T Qkx(t) is equal to —x(t) T x(t), 
then the minimum rate must be equal or smaller than 
the convex combination, and lemma is proved. □ 


This Lemma ensures that a generator Aj that makes 
the Lyapunov function (12611 strictly decrease always ex¬ 
ists, and it provides also indications on which system 
allows a fastest decrease of the Lyapunov function. 


B. Steepest-descent switching 
with fixed switching interval 

As suggested by the Lemma of the previous section, a 
switching law that attains the fastest possible decrease 
rate of the Lyapunov function that is associated to the 
stabilizing convex combination A c can be simply ob¬ 
tained by selecting, at each time, the generator that en¬ 
sures the steepest descent. Namely, the current value of 
the switching signal j(t) is chosen as follows: 

j(t) = argmi n{xjQ k x t }. 

kdV 

If two or more generators correspond to the same de¬ 
crease rate, we can pick either one. 

In order to (locally) maximize the convergence rate, 
the switching signal should switch as soon as any gener¬ 
ator allows a faster decrease of the Lyapunov function - 
that is, one should compute the steepest-descent evolu¬ 
tion at all times. 

This approach has two issues: the first one is that it 
may not be possible to implement such switching law 
- it requires to either to be able to determine the full 
(continuous-time) switching law off-line, or be able to 
have some mapping of the state space in areas where a 
generator is optimal, and perform the switching when 
the current state, obtained by integration of the imple¬ 
mented dynamics, reaches a border between the areas. 
The second issue is that we do not know how to guaran¬ 
tee that a finite number of switches occurs in any finite 
time interval. Similarly, we do not know how to guar¬ 
antee convergence if the switching time is bounded for 
below. In cases where the computational task of deter¬ 
mining the optimal switching law is not viable, or we 
want to prevent an infinite number of switches from oc¬ 
curring in a finite interval, suboptimal solutions can be 
devised. 

The simplest way to bypass the problems described 
above is to impose a minimal switching time AT > 0, 
and compute the steepest-descent generator at the be¬ 
ginning of each interval. On the one hand, this of course 
prevents infinite switching in any finite-length interval. 
On the other hand, by Lennna [1] along with the continu¬ 
ity of the evolution and its derivatives (apart from the 
switching instants), we know that there exists a small 
enough AT > 0 such that the Lyapunov derivative re¬ 
mains strictly negative for all times. This implies con¬ 
vergence of the steepest-descent method for some mini¬ 
mal switching interval. However, in some cases the latter 
may be to small or hard to compute. Another suboptimal 
strategy is described in the following section. 


C. Sub-optimal switching 
with minimal guaranteed descent 

A suboptimal solution with a guaranteed rate of con¬ 
vergence can be devised as follows [29|. Let us define a set 
of arbitrary real numbers rj £ (0,1], each of which is as¬ 
sociated to a generator Aj. These are tuning parameters 
for the algorithm and can be taken smaller than one in 
order to avoid chattering and overly fast switching laws, 
at the price of slower convergence. They correspond to 
the minimal accepted descent rates for V. If all rj are 
taken to be one, the convergence is at least as fast as the 
one associated to the non-switching dynamics A c . 

The algorithm is initialized as the steepest descent one: 
at the first instant to a j is selected such that: 

.j(f 0 ) = a.rgunn{xf 0 Q 1 x to ,..., xf 0 Q m Xt 0 }- 

The next switching instant is then chosen to be 

ti = inf{f > t 0 : x T (t)Q j{to) x{t ) > -r j{to) x T {t)x{t)}, 

and the following active subsystem is that associated to 

j{ti) = argmin{a:(ti ) T Qix{t\),... ,x{ti) T Q m x{ti)}. 
j£V 

By Lemma [TJ we are guaranteed that such a generator 
ensures that the Lyapunov function strictly decreases, 
and does so faster than it would using the stabilizing 
convex combination. The sequences of switching times 
and active systems can be thus recursively defined: 

tk+i = hif{f > t k : x T {t)Q j{ t k )x{t) > -rj( tk )X T (t)x(t)}, 

(27) 

and 

j{tk+ 1 ) = argmin{x(tfc+i) T Qi^(ifc+i), • ■ ■ (28) 

p£V 

. . . , x(tk-\-l ) Qm^fffc+l )} 1 

Notice that in order to compute the switching law, we 
need the current state (either by off-line calculations, or 
real-time simulation in the spirit of model-based feed¬ 
back El)- In this case, we only solve a minimum problem 
when the descent becomes slower than the active thresh¬ 
old rj. 

Given such a switching law, we need to check whether 
the switching system is well-posed, that is, if the set of 
jump times is finite for any finite interval. 

Theorem 2 Under the above switching law, the switch¬ 
ing system is well-posed and asymptotically stable. 

The proof, that follows known results for classical sys¬ 
tems [29}, is provided in Appendix lEil for completeness. 

D. Robustness with respect to initialization 

In principle, one needs to know the initial state of the 
quantum system of interest in order to apply state-based 
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switching techniques. An accurate estimation of such 
state appears to be, in general, necessary for the success¬ 
ful stabilization of the target state. 

Example: Let’s suppose, for example, that the target 
state is 


P = 


1 0 0 
0 0 0 
0 0 0 


Since supp(po) Q supp(po) by hypothesis, for any such 
density operator po, there exists a A > 0 such that, 


Po < Apo- (32) 

In particular, one can always take A to be the inverse 
of the smallest eigenvalue of po- Since the evolution is 
positivity preserving, it preserves matrix ordering and 
we can then conclude that 


and that we dispose of two Lindblad generators, respec¬ 
tively associated to the following noise operators: 


lim £ t {po) < A p. (33) 

£—>+oo 


o 

T- 1 

o 


0 0 o' 

0 0 0 

l 2 — 

0 0 1 

0 0 0 


0 0 0 


It is easy to see that they satisfy Assunption 1: any non- 
extremal convex combination is in fact stabilizing - this 
can be verified easily with the methods of mm- if the 
estimated initial state is 


Po 


0 0 0 
0 1 0 
0 0 0 


then the state-based algorithms will try to use only the 
first generator to reach the target faster - in fact the sec¬ 
ond generator would leave the estimated state invariant. 
But if the actual state is instead 


Po 


0 0 0 
0 0 0 
0 0 1 


then the first generator alone does not stabilize the tar¬ 
get, as po is an invariant state for it. □ 


This example proves that if the switching rule is de¬ 
signed for the wrong initial state, then the algorithm may 
fail to provide a stabilizing law. However, we next present 
a result that allows us to design state-based strategies 
that are robust with respect to variations of the actual 
initial state, as long as the estimated one is a full rank 
state. 


Proposition 5 If a family of TPCP maps {£t}t^o 
such that 


lim £ t (p 0 ) = p, (29) 

£—>•+00 


then for any po such that supp(po) C supp(po) we have 
that 


lim Tr(nf t (p 0 )) = 1, (30) 

t— >+oo 

where n is the orthogonal projector on supp(p). 

Proof. As £t is linear and the limit (l29l) exists, then for 
any scalar A we have 

lim £ t (Xp 0 ) = Xp. (31) 

t— >+oo 


Using m and the fact that the evolution is trace¬ 
preserving we have: 


lim Tr(n£ t (p 0 )) = lim Tr(£ t (p 0 )) = 1- 

£—>+oo t— >+oo 


□ 

Assume now p to be a pure state. Then n = p is a one- 
dimensional projector, and the only density operator for 
which A p is an upper bound is p itself. Furthermore, 
the condition supp(po) Q supp(po) is satisfied if po is 
assumed to be full rank. We thus have the following: 

Corollary 2 If a family of TPCP maps {£t}t^o such 
that 


lim £ t (p 0 )=p, (34) 

t— >+oo 

for a full-rank initial state po and a pure target state p, 
then p is GAS, i.e. it is asymptotically reached indepen¬ 
dently of the initial state po ■ 

It is worth noting that the same results still holds if we 
relax CP to simple positivity, or if we consider discrete¬ 
time evolutions [10| . 

These stability results immediately imply that if a 
state-based switching law is designed for a full-rank ini¬ 
tial state po in order to prepare a pure target state p , 
then the corresponding evolution makes p GAS, i.e. it is 
asymptotically reached independently of the initial state 
Po • Analogous reasoning can be obviously applied to sub¬ 
space preparation. 

In designing the control law we can always consider 
the completely mixed initial state, as one would do in 
the absence of initial information, or, if we have a good 
estimate po which is not full rank, we can consider a 
perturbed version po = (l — e)po + eI/d, with e arbitrarily 
small, in order to ensure robustness. 

We could expect that choosing a state-based switching 
rule, without exploiting the information on the initial 
state, would lead to a stabilizing law characterized by a 
convergence rate comparable with a time-based switch¬ 
ing rule, for which the initial state is a priori supposed 
to be unknown. Nevertheless, in the numerical examples 
that follow, we have that it still converges to the target 
with better performances than those obtained by imple¬ 
menting a time-based switching law. 
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V. EXAMPLES 

In the following we compare the performances of dif¬ 
ferent Markov evolutions that prepare an entangled pure 
state of interest, in terms of both a natural Lyapunov 
pseudo-distance and the Euclidean distance in coherence- 
vector representation (equivalent to the Hilbert-Schmidt 
norm). The evolutions we compare for both examples are 
four: 

1. A time-independent semigroup generator, associ¬ 
ated to the stabilizing convex combination of As¬ 
sumption [TJ for which the target p is known to be 
the unique equilibrium (denoted in the figures as 
no-switch)', 

2. A cyclic, time-based switched evolution, such that 
a convex combination of the alternating generators 
is equivalent to time-independent generator above 
(i.e. they satisfy Assumption 1). For both ex¬ 
amples, a convex combination with uniform coeffi¬ 
cients is sufficient, and used. A minimal switching 
time AT is chosen, so that the evolution is stabiliz¬ 
ing. The Lyapunov function (l26l) associated to the 
stabilizing convex combination will be called the 
natural Lyapunov function', 

3. A state-based switching evolution, where the gener¬ 
ator is chosen at AT intervals in order to guarantee 
the fastest decrease of the natural Lyapunov func¬ 
tion, as described in Section HV 61 depending on a 
full-rank estimate of the initial state. The same 
evolution is also applied to a different initial state, 
in order to verify the robustness of the approach, 
and the potential deterioration in performances; 

4. A state-based switching evolution, where the gen¬ 
erator is chosen as described in Section II V Cl de¬ 
pending on a full-rank estimate of the initial state. 
The same evolution is also applied to a different 
initial state, in order to illustrate the robustness 
of the approach, and the potential deterioration in 
performances; 

Choosing a pure target allows us to compare the per¬ 
formance of state-based switching strategies with correct 
initializations and faulty ones, since the results of Sec¬ 
tion [IVD] guarantee convergence. Using a distance that 
is not tailored to the problem (the Euclidean one) help us 
to illustrate how the convergence in Lyapunov distance 
does not hides undesired behaviors. For all simulations, 
we use an integration step (the time unit for the X axis of 
the plots) of t = 0.02, and a minimal switching interval 
AT = 3 1. 


A. Bell states 

We first consider a two-qubit system defined on Hab = 
Ha ®Hb ~ span{|0), |1)}® 2 . Our aim is to prepare the 


maximally entangled Bell state, 

Pab = ^ (100) + |11»«00| + (11|), (35) 

by switching between Lindblad dynamics and to compare 
time-based and state-based switching rules convergence 
rates. Denote as usual the Pauli matrices as: 



'0 

1' 


' 0 i 


'1 0 ' 

&X — 

1 

0 

j a y ~ 

-i 0 

, <Tz = 

0 -1 


with the scaled identity ao = ^=J. 

The first generator we consider is associated to the 
Hamiltonian H = a y <g>cro+<Jo ®tr y . The second generator 
is determined by the Lindblad operator L = a z <S> cr 0 — 
i(a y <g> cr x ). The target state pab is a fixed state for 
both the generators and it can be proved that applying 
them simultaneously, with any nonzero weight, makes 
Pab asymptotically stable [tJ . 

Let us denote L i and L 2 the superoperator form of the 
Lindblad generators derived by H and L respectively, 
according to the coherence-vector representation, as in 
III Bl After applying the change of basis that makes the 
generators of the reduced dynamics linear, it is easy to 
numerically check that A c = + \A 2 is actually a 

Hurwitz convex combination. This combination satisfies 
Assumption 1 and allows us to derive time-based and 
state-based switching rules that prepare the target state 
as described in the introduction to this section. 

In Figure [Hand [2] we compare the convergence features 
in terms of the natural Lyapunov function, and in terms 
of the coherence-vector euclidean distance to the target, 
respectively. The actual initial state is p = |0)(0| (g) |0)(0|, 
while the estimated initial state is chosen to be the com¬ 
pletely mixed one. The estimated state is used to com¬ 
pute the state-based switching strategies, whose esti¬ 
mated behavior is depicted with dashed lines. All the 
solid lines represent the evolution of the distances calcu¬ 
lated for the actual states. We set the parameters for the 
suboptimal state-based switching to be r± = T 2 = 1. 

In Fig[T] we find, as expected, that the time-based 
switching is the slowest to converge, and at close inspec¬ 
tion it even exhibits local increase in the Lyapunov func¬ 
tion. In fact, time-based switching is an approximation 
of the no-switching strategy, and during some switching 
steps the state may also tend to get farther from the tar¬ 
get - it is only after a full cycle that we are guaranteed, 
for sufficiently fast switching, that we the distance to the 
target has decreased. In our simulations, the bad perfor¬ 
mance is caused by the fairly large cycle time, making the 
cyclic switching a rough approximation of the reference 
convex combination. In this particular case, the effect is 
likely strengthened by the presence of a purely Hamilto¬ 
nian generator among the switching ones, which does not 
induce dissipative, and hence contractive, effects. 

The convex-combination, Hurwitz generator displays a 
smooth converging behavior, with the distance from the 
target state with the state-based rule is monotonically 
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Time (Integration Units) 


appears to be faster for the actual state. 


B. GHZ states 


We here present an example of preparation of a three- 
qubit GHZ state, 

Pghz = 1(|000> + |111»«000| + (111|), (37) 

with the switching techniques we introduced above. The 
stabilizing, non switching generator we consider is one 
that has been proposed in [38|], and can be implemented 
via an Hamiltonian 


H = a « - ai 2 ) 


rO) 


FIG. 1: (Color online) Evolution the Lyapunov function with 
different switching rules in the two-qubit case. 



Time (integration Units) 

FIG. 2: (Color online) Evolution the Euclidean distance with 
different switching rules in the two-qubit case. 


decreasing. The state-based strategies show, for both 
the estimated and the actual state behavior, faster con¬ 
vergence. The evolution of the actual state converges to 
the target state but it is not monotone, as it is clear by 
looking at the peaks around 50 time units. Nonetheless, 
it remains much more effective than the evolution asso¬ 
ciated to convex combination of Assumption [l] where no 
switching is performed and all generators are activated at 
each time scaled by the corresponding weights, for both 
state-based strategies. The “steepest descent” approach 
is, as expected, overall slightly faster than the suboptimal 
one with guaranteed minimal descent. 

The qualitative evaluation of the plots of Fig [2] leads 
to the same conclusions, confirming the advantage of the 
state-based strategies, even in presence of a wrong initial¬ 
ization of the algorithm. In fact, the initial convergence 


where a i*'* = I <g) ■ ■ ■ <g> a x <g) • • ■ <g> /, and a x acts on the 
*-th qubit, along with two noise operators, 

L x = (|00)(01| + |11)(10|)®/, 

i 2 = J®(|00)(01|+i|ll)(10|). 

The simultaneous action of the three components, with 
arbitrary positive weights, leads to the asymptotic stabi¬ 
lization of pghz- A stabilizing time-based switching rule 
can be therefore easily obtained with a fast periodical 
switching between them. We then consider switching be¬ 
tween £o, 1 , 2 , corresponding to the superoperator action 
of H, Lx,L 27 respectively. 

The actual initial state is p = |0)(0| ® |0)(0| ® |0)(0|, 
while the estimated initial state is chosen to be the com¬ 
pletely mixed one. The estimated state is used to com¬ 
pute the state-based switching strategies, whose esti¬ 
mated behavior is depicted with dashed lines. All the 
solid lines represent the evolution of the distances calcu¬ 
lated for the actual states. We set the parameters for the 
suboptimal state-based switching to be rr = r 2 = = 1. 

In Fig|3] we find, as expected, that the time-based 
switching is the slowest to converge in terms of the Lya¬ 
punov function, and close inspection reveals that the cy¬ 
cling indices some local fluctuations. This is because 
of the fairly large cycle time, making the time-based 
switching a rough approximation of the reference con¬ 
vex combination. The effect may also be partially due 
to the presence of a purely Hamiltonian generator that 
does not induce dissipative contraction effects. The non¬ 
switching generator associated to the stabilizing convex 
combination displays a smooth, monotonically decreas¬ 
ing behavior, as expected. The state-based strategies 
show, for both the estimated and the actual state be¬ 
havior, the fastest overall convergence. The evolution 
of the actual state converges to the target state but the 
Lyapunov pseudo-distance does not start decreasing im¬ 
mediately, a clear sign that the optimal generator for 
the completely mixed state is not effective for the actual 
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FIG. 3: (Color online) Evolution the Lyapunov function with 
different switching rules for the three-qubit case. 



FIG. 4: (Color online) Evolution the Euclidean distance with 
different switching rules for the three-qubit case. 


initial state. Nonetheless, after another generator is se¬ 
lected, state-based strategies become more effective than 
the evolution associated to the convex combination of As¬ 
sumption [0 without switching. The “steepest descent” 
approach is overall slightly faster than the suboptimal 
one with guaranteed minimal descent. 

A careful scrutiny of the plots of Fig 01 for the undis¬ 
torted euclidean distance in coherence vector space, es¬ 
sentially leads to the same conclusions, confirming the 
advantage of the state-based strategies, even in presence 
of an incorrect initialization of the algorithm. 


VI. DISCUSSION AND PERSPECTIVE 

In this paper we have presented two approaches for de¬ 
signing switching quantum Markov evolutions that make 
a desired state GAS. For algorithms obtained by using 
either time-based or state-based designs, the underlying 
assumption is the existence of a convex combination of 
the switchable generators that would make the target 
GAS. The first approach we propose constructs time- 
based, cyclic switching laws: its advantages are the sim¬ 
plicity, and the intrinsic robustness with respect to the 
initial state. No initial state estimate, optimization or on¬ 
line computation is needed. However, being essentially a 
way to approximate the stabilizing convex combination, 
the method is quite sensitive to the minimum allowed 
switching time. Numerical simulations, as exemplified 
by those presented in this paper, show that this strategy 
is in general the worst performer of the group. 

State-based strategies, while they entail some com¬ 
putational overhead, guarantee to outperform the non¬ 
switching evolution when the initial state is accurately 
known, and numerical simulations show that they retain 
this advantage even with a generic initialization (e.g. the 
completely mixed state we used). Robustness with re¬ 
spect to the initial state can be guaranteed if the target 
is a pure state (or a subspace) by arbitrarily small pertur¬ 
bation of the initial estimate. The state-based strategies 
thus offer a valid alternative to time-invariant methods 
for dissipative preparation of states, especially when im¬ 
plementation of complex dynamics in an experimental 
setting is challenging, and further motivates their study. 

The structure of the algorithm suggests that bet¬ 
ter results, in terms of robustness, may be attained 
by combining optimal state-based switching with state- 
estimation methods. The latter could in principle be 
implemented via continuous measurements and filtering 
equations, when these guarantee asymptotic convergence 
of the estimated state to the actual one [Tsl . j32], [§9}. It 
is known that, for pure states and subspaces with time- 
independent stochastic dynamics, conver gen ce in expec¬ 
tation implies convergence in probability [32|. An exten¬ 
sion of these results would help proving the effectiveness 
of the filtering-based switching. This combination would 
lead to a control strategy that closely resembles a classi¬ 
cal one, where state-based switching is typically decided 
via evaluation of the current state through a feedback 
loop. 

Other potential developments of these results may in¬ 
clude a comparison with the performance of optimal con¬ 
trol methods for state preparation and integration with 
scalable strategies for the dissipative preparation of en¬ 
tangled states on large networks [3(| • 
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Finally, by denoting: 

«i = 3?? Il^( 0 >*)ll> 

0<t<£ 

which is always attainable because 0(0, t) is continuous 
in t, we get 


Appendix A: Proof of Proposition [l] 


||0(s 2 ,si)|| < n\n.e- x{l2 - h)e < ^ Ke - A(s2 - Sl) . 


Proof. Consider a periodic switching with period e. 
During each period, each system Aj acts for a time in¬ 
terval proportional to its coefficient aj (the order is not 
important, as we shall see). The transition matrix over a 
period can be expressed as the exponential of an effective 
generator A: 

exp(Ae) = exp (a m A m e)... exp(ai A\e). 


We can then separate the contribution of the average 
generator A c with respect to the part depending on e as 
A := A c + T c e and T c is bounded and contains all the 
higher-order terms of a Taylor expansion of A around 
e = 0. The eigenvalues of a matrix continuously depend 
on its entries. As Y c is bounded, if e —> 0, then the 
eigenvalues of A approach those of A c and assuming that 
A c is Hurwitz, there exists a e > 0 for which A is Hurwitz 
too. Fixed such an e, a periodic switching path can be 
defined e.g. as follows: 


1 

2 

cr(t) = < . 


if mod (f, e) £ [0, are) 
if mod (f, e) £ [ai£, a\£ + a 2 e) 

(Al) 


m if mod (t,e) £ [(J]’’0 1 1 a j )e,e) 


The transition matrix is exponentially convergent and 
the switching system is then stabilized. □ 


Appendix B: Proof of Theorem [2] 

Proof. We want to show now that such a switching law 
is acceptable, that is actually tk+i > t k - Let 6 be any 
real number greater than 1. We first consider the case: 

lk(i)|| < Q\\xk+i\\ Vf € [ffc,f fc+ i]. (Bl) 


and define 

g{t) = x(t) T (Qj + I)x(t) t £ [tk, tk+ 1 ], 

where j = cr{t k +). On the one hand, thanks to LemmaU 
we know that 

x{t k { QjX{tk) — x(tk) x{t^ 7 
and hence we are guaranteed that 

g(t k ) < 0. (B2) 

On the other hand, from (1271) ,we have 


Setting si < S 2 , let us define 

0( S2 , Sl ) ;= e A p(s2-ip) e a P -iAp_ie _ _ _ e A fc _i(t fc -«i) ^ 


gitk+i) > (1 - rj)xl +1 x k+ i. (B3) 

Deriving the latter we get 


the transition matrix from the state at instant si £ 
(t k -i,t k ) to that at instant S 2 £ (t p ,t p + 1 ) according 
to a given switching law. For any non-negative integers 
h < h the evoution covers a finite number of cycles, so 
that 

cj)(l2£, he) = 

As A is Hurwitz, there exist positive numbers k and A 
such that 

\\<j>(he,lxe)\\ < Ke~ x{l2 ^ h)s . 

For any si < S 2 , let ii and I 2 satisfy 

he < si < (l 1 + l)e, {h - l)e < s 2 < he. 


f t = x(t) T (Aj (Qj + I) + (Qj + I)Aj)x(t). 

By denoting 

Vj := II Aj {Qj + I) + {Qj + I)Aj\\, 
and using (EH), we have 

|^| < d 2 r]jX k+1 Xk+i Vf £ [t k ,tk + !]■ (B4) 

According to (IB2I) and (IB3I) 

g{t k+1 ) ~ g{t k ) > (1 - r 3 )xl +1 x k+l 
tk- f-1 t k t k -\-1 t k 


Then we have: 


Remembering then (IB4I) . 


|| 0 (s 2 ,si)|| < \\<t>{he, si)11 • ||0(Z 2 e,Zie)|| • ||0(s 2 , he)\\ 
< «:e _A(i2£_ ' l£) ||0(O, si - Zie)|| • ||0(O, l 2 e ~ s 2 )||. 


( l ~ r j) X k+\X k +i 
tk -\-1 tk 


< 9 2 VjX k+ ix k +i, 
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and, consequently, 


Finally, gathering both the cases, we can say that 


tk+i ~ tk > 


(! ~ r j) 
d 1 2 3 4 5 6 7 8 9 10 11 rjj 


This shows that the switching instants are have mini¬ 
mal distance between them, and there cannot be infinite 
switching times in any finite time interval. 

We are then left with the other case, when m does 
not hold, and thus: 


3 1* G [t k ,t k+1 ) : ||x(t*)|| > 0||a;fc + i||. 


t k -(-1 — t k > sup min 


1 - Ti 


Hd) \ 


0 2 Vj '\\AjW)' 


and that the switching signal is valid as the difference is 
always positive. 

Choosing then V(x) = x T Px as Lyapunov function we 
get 


As the system dynamics in this time interval is described 
by Aj. then 


x(t*) = exp (A p (t* - t k+ i))x k+1 . 
From the latter, and remembering that 


dV 

dt 


x T (t)Q a ( t )x(t) < r rT ( t ' ) x T (t)x{t) < - rx T (t)x(t ), 


where 


II exp(Aj(t* - 4+i))|| < exp(||Aj||(t fc+1 -t*)), 
it follows that a minimal spacing is also guaranteed: 

ln(0) 


r := min{n,... ,r m }, 


tk+l ~~ tk > tk+1 — t* > 


and the theorem is proved. 


□ 
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